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^ Abstract 

O A difficulty with previous treatments of the gravitational self-force is that an explicit formula 

<N 

for the force is available only in a particular gauge (Lorcnz gauge), where the force in other gauges 

^ must be found through a transformation law once the Lorenz gauge force is known. For a class of 

^ gauges satisfying a "parity condition" ensuring that the Hamiltonian center of mass of the particle 

,— I is well-defined, I show that the gravitational self-force is always given by the angle-average of the 

O 

bare gravitational force. To derive this result I replace the computational strategy of previous work 

with a new approach, wherein the form of the force is first fixed up to a gauge-invariant piece by 

^ simple manipulations, and then that piece is determined by working in a gauge designed specifically 

> 

to simplify the computation. This offers significant computational savings over the Lorenz gauge, 

cn 

\^ since the Hadamard expansion is avoided entirely and the metric perturbation takes a very simple 

in 

form. I also show that the rest mass of the particle does not evolve due to first-order self-force 

o 

effects. Finally, I consider the "mode sum rcgularization" scheme for computing the self-force 

^ in black hole background spacetimes, and use the angle-average form of the force to show that 

^ the same mode-by-mode subtraction may be performed in all parity-regular gauges. It appears 

H 

plausible that suitably modified versions of the Regge- Wheeler and radiation gauges (convenient 
to Schwarzschild and Kerr, respectively) are in this class. 
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The leading-order deviation from geodesic motion proportional to the mass of a body — 
interpreted as the force due to the body's own gravitational field — is known as the grav- 
itational self-force. A recurring source of difficulty in both the theoretical treatment and 
practical computation of the self-force has been the choice of gauge in which the metric per- 
turbation and force are expressed. In particular, early treatments P-H] required a specific 
gauge choice — Lorenz gauge — even to define the perturbed trajectory (via a point particle 
hypothesis coupled with "Lorenz gauge relaxation" |2l[5] to allow non-geodesic motion), and 
a proposed extension of the results to other gauges [Hj restricted to gauge vectors that are 
continuous at the particle, even though the metric perturbation is singular. At a theoret- 
ical level, this elevates a particular gauge to fundamental status, to the exclusion of other 
gauges that seem equally nice, such as any gauge where the point particle 1/r singularity 
corresponds to linearized Schwarzschild in Cartesian Schwarzschild coordinates, as opposed 
to the Cartesian isotropic coordinates that correspond to Lorenz gauge. ^ And at a practical 
level, one has excluded the standard gauges of black hole perturbation theory [B] . 

Previous work [S] (hereafter paper I) eliminated the fundamental status of the Lorenz 
gauge by giving a definition of perturbed motion holding for any gauge where the particle 
is represented by a 1/r singularity. However, in this work we still relied on the Lorenz 
gauge for our computations and, more importantly, for the expression of the final result, as 
a formula holding in the Lorenz gauge together with a generalized transformation law. At a 
theoretical level, we still have a preferred role for the Lorenz gauge; and at a practical level, 
the results suggest that the computation of a self-force in an alternative gauge must always 
proceed through Lorenz gauge, eliminating much of the appeal of working in alternative 
gauges in the first place. 

In this paper I will identify a class of gauges based on the requirement that the center 
of mass as defined by Regge and Teitelboim [7] is well-defined (in the "near zone"), and 
show that the force in any such gauge is given by the angle-average of the bare force in that 
gauge. To derive this result I adopt the assumptions of paper I but take a new computational 
approach, wherein the form of the force in any gauge is fixed up to a gauge-invariant piece 
by simple manipulations, and then that piece is determined by working in a gauge chosen 

^ The (discontinuous) gauge vector that changes the singularity from isotropic to Schwarzschild type is 
given by — — x^/r in Fermi normal coordinates. 
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specifically to make the computation as simple as possible. This approach avoids much 
of the computational complexity of previous work (eliminating the Hadamard expansion 
entirely and significantly reducing the calculation needed thereafter), while organizing the 
computation so that the final equation automatically takes a gauge-independent form. 

The precise results are as follows. We define the (lowest-order) mass M, spin Sab, and 
center of mass deviation of the particle as tensors on a timelike worldline 7 (four velocity 
u"') in a vacuum background metric Qab-"^ Then we find that that 7 is a geodesic, that the 
mass and spin are constant (parallel propagated), and that the deviation Z°- satisfies 

u'VbiW^V^Z-) = ^ lim I F-dQ + R.Ju'Z^u" + ^R^JS'^u\ (1) 

where F"" is the "bare gravitational force", 

F- = - {g-^ + u'^u') (Vdhc - l^bKd^ u'^u\ (2) 

and h^j^y is the metric perturbation of a point particle, 

Gi;) [h] = SttM I 6^{X, z{T))UaUbdT, (3) 

J y 

which must be expressed in a "parity-regular" gauge, where the singular part of the spatial 
metric is even-parity on the sphere, Cij{t, —n) = Cij{t, n) in equation (|6]). The bare force F°' 
is familiar from the perturbed geodesic equation, and here is defined only off of 7 (and only 
locally, where u"" is extended off 7 by parallel transport along spacelike geodesies orthogonal 
to m") since the metric perturbation is divergent. The integral in ([T]) is defined by using 
the exponential map based on 7 to associate a fiat metric, in terms of which the integral 
is over a fixed 2-sphere of spatial distance r in the hyperplane orthogonal to u", with ra" 
its unit normal and dVt its area element, and the integration is done component-wise under 
the exponential map. I also define the perturbed mass of the particle and show that it is 
constant in time. 

The first term on the right hand side of equation ([T]) is proportional to the metric per- 
turbation; I therefore refer to this term as the gravitational self-force. We see that the force 
in any parity-regular gauge is given by the angle-average of the bare force in that gauge, 
so that the self-force may be viewed as the net gravitational force on the particle. If the 

^ The spin is antisymmetric and satisfies u°'Sab = 0. The deviation satisfies Z'^Ua ~ 0. 
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Lorenz gauge (which is parity-regular) is adopted and the Hadamard form for the metric 
perturbation is computed (choosing the retarded solution with no incoming radiation) and 
plugged in, then this term reduces to the standard "tail integral" expression for the self-force 
(e.g., [3]). The second term corresponds to the geodesic deviation equation and reflects the 
particle's desire to move on a new geodesic once it has been displaced from the original. 
The third term is the Papapetrou spin force. If the parity condition is violated, equation ([T]) 
does not hold, and the equation of motion takes a complicated form involving an explicit 



gauge transformation to a parity-regular gauge, equation (B2). 

A practical technique for computing self-forces in black hole background spacetimes is 
known as "mode sum regularization" (|8j and many other references). In this approach one 
numerically solves for the spherical harmonic modes of the metric perturbation (the sum over 
which diverges at the particle), and performs a mode-by-mode subtraction that regularizes 
the sum in such a way that the correct self-force is computed. Extensive work has determined 
the form of the subtraction in the Lorenz gauge, which automatically holds for gauges that 
are smoothly related. Taking advantage of a connection between mode decompositions and 
averaging (and hence self- force), I show that the same subtraction may in fact be performed 
in all parity-regular gauges. It appears plausible that suitably modified Regge- Wheeler and 
radiation gauges (convenient to perturbations of Schwarzschild and Kerr, respectively) will 
be in this class, with especially strong evidence in the case of radiation gauge (see discussion 
at the end of section |V]). If so, then the results of this paper would provide a theoretical and 
practical foundation for the computation of self-force effects in gauges convenient to black 
hole spacetimes. 

I use the conventions of Wald [11]. Greek indices label tensor components, while early- 
alphabet Latin indices a,b, . . . are abstract indices. When coordinates t, are used, the time 
and space components are denoted by and mid-alphabet Latin indices respectively. 



II. REVIEW OF FORMALISM 



The central idea of paper I is to introduce a mathematically precise formulation of the 
notion of the "near-zone" of a body, and to use the requirement of a sensible near-zone 
to demand a sensible perturbation family. Given a family of metrics Qabi^) in coordinates 
(t,x*), we define a scaled metric (jab = ^~'^gab and scaled coordinates t = A"^(t — to) and 
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X* = \~^x\ Denoting the scaled metric in scaled coordinates by (jfip, consider the A — )■ 
limit, gj^p = giip\\=o- This limit effectively "zooms in" on the spacetime point {t = to, = 0), 
and will recover the near zone of a body if the one-parameter-family contains a body whose 
radius and mass shrink down linearly with A to the worldline = (denoted by 7). By 
demanding its existence (and associated conditions), we automatically consider bodies of 
small size and mass. Note that the components of the original and scaled metrics are 
related simply by "plugging in the new coordinates," 

QfipiX; to; t, x') = g^^{X; t = to + At, x' = Ax*). (4) 

This equation relates components of the scaled metric in scaled coordinates, g/ip, to corre- 
sponding components of the original metric in the original coordinates, g^j^^,. 

One may perform perturbation theory in either the original ("far-zone") picture or the 
scaled ("near-zone") picture. The far-zone background and perturbations will be denoted 
by gjj]) = {l/n\){dx)"'gfj^„\x=o (I also write hf^i, = gjxJ), while the near-zone background 
and perturbations will be denoted by gj^J = {l/n\){dx)"'giip\x=o. Our assumptions together 
with the choice of Fermi coordinates (e.g., [4J) about 7 for the far-zone background metric 
constrain the far-zone quantities to take the form 

= %u + B,,,,it)x^x^ + 0{r') (5) 
V = = + D,,it, n) + rE,,{t, n) + 0{r') (6) 

9t} = + + K,At, n) + 0{r\ (7) 

where we have defined r = SijX^x^ and ra* = x^/r. Here -B^i^j is related to the Riemann 
curvature on the background worldline x* = (see e.g. j4j for the exact expression), whereas 
C^^, . . . , Kf^i, are unspecified smooth functions on M x 5*^ (denoted by arguments {t, n)). The 
lack of a linear term in equation (|5| is a consequence of 7 being geodesic, as shown in paper 
I. We also showed that the metric perturbation has effective point particle source, equation 
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(|3|. In light of equation ^ the near-zone perturbation series takes the form 

1 
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where an overdot denotes a t-derivative, and the order symbols 0(l/f") refer to t- 
independent functions. The dependence on {to,n) is suppressed in equations ^ and (10) 
for readability. Note that the lack of growing-in-f terms in equation ^ is inherited from 
the absence of a linear term in ([s]), which is a consequence of 7 being geodesic (and choosing 
Fermi coordinates). The (stationary and asymptotically fiat) near-zone background metric 
g'pp represents the local state of the body at time to, and as such its properties should char- 
acterize those of the body. In particular, the multipole moments of g^^p should correspond 
to those of the body, at lowest non-trivial perturbative order. To the (far- zone) perturbative 
order considered in this paper, only the monopole and dipole moments can play a role. We 
therefore define^ 

(11) 
(12) 

(13) 

and refer to M, D\ and Sij as the (lowest-order) mass, mass dipole, and spin (current dipole) 
of the body."^ In equations (ll)-(n!3|, the bars on the indices of coordinate components of the 
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1 n^g^^Jf^dn 
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f gSln'f^dn 




3 

" 8^ 


lim 

r~^oo J 


fg^^^n^f^dn, 



near-zone metric have been dropped, and the integrals are taken on fixed coordinate two- 
spheres with respect to the "flat" volume element dQ. It is well known (and easily checked) 
that the mass dipole may be set to zero by translating the coordinates, x* — )■ a;* — Z*, with 

Z\to)=D\to)/M{to). (14) 



Equation ( 13 1 for the spin holds only in coordinates where C^o vanishes. A formula for the spin holding 



in general coordinates would take a more complicated form. 

The mass dipole may also be thought of as the time-space component of the spin tensor, S''^'. However, 
we will work with a spin tensor that is orthogonal to 7, 5°* = 0, defining a separate mass dipole. 
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This gives the new coordinates the interpretation of being mass centered, so that rep- 
resents the center of mass position of the body in the original coordinates. Since the 
above translation corresponds to a first-order gauge transformation in the far zone (recall 
X* = xYA), we identify with the first-order deviation of the center of mass position from 
its background position 7. While the mass, spin, and deviation are defined in the near zone 
"at spatial infinity" (as a function of the coordinate time to along 7 at which the near-zone 
limit is taken), they may equally well be viewed as tensors defined in the far- zone "at the 
spatial origin", i.e., as tensor fields on 7. It was shown in paper I that the mass and spin 
do not evolve with time to (i-e., that they are parallel propagated tensor fields on 7), while 
a Lorenz-gauge equation of motion (together with a gauge transformation law) was derived 
for the deviation vector. 

The analysis of this paper will require transformation properties of not only under 
ordinary translations, but also under transformations of the form bx^ = a*(n) +0(l/f), i.e., 
under angle-dependent translations, or supertranslations. Using the well-known fact that 
goQ = — 1 -|- 2M/f + 0(l/f^) for all vacuum solutions of the form ([s]), we have by direct 
computation that 6goo = 2Ma^ni/f'^ + 0{l/f^), so that transforms as 



Since the angle-average of a vector picks out an £ = 1, electric parity part, we may restrict 
consideration to a* of the form a* = B^njn^ + for constants -B* and C*. The associated 
change in is then simply + C". Thus changes under supertranslations with 5* 7^ 0, 
in addition to its change under ordinary translations a* = C\ 

III. COORDINATE CHOICES AND PERTURBED MASS 

The mass, spin, and deviation of the body were defined at lowest-order in the near 
zone, so that perturbative corrections to these quantities would naturally be defined at first 
and higher orders in the near-zone. However, while the background metric is stationary 
and asymptotically flat (so that its multipole moments are well-defined), the n'th-order 
perturbation may have growth in n combined powers of t and f (c.f. paper I and equations 
(|8]|9])). It turns out, however, that at first order the situation is more under control. As 
already noted, the choice of Fermi coordinates in the far-zone background — together with the 




(15) 
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fact that 7 is geodesic — eliminates growing-in-f terms from the first near-zone perturbation, 
so that the perturbation is asymptotically fiat. Furthermore, as I now show, the constancy 
of the mass and spin similarly allows one to eliminate all i/f and t/f^ terms, so that the 
perturbation is stationary to 0(l/f^). To see this, recall that g'^^^ is a stationary and 
asymptotically fiat metric, so that, introducing 5^,^ = diag(l, 1, 1, 1) and = (1,0,0,0), it 
may (at each to) be put in the standard form |13j . 

^ = ... + + |;(3,„„ - UK) - + O (i) , (16) 

where SiQ = for all to (implying = 0, so that the coordinates are mass centered and 
"track" the motion of the body), while a priori M and Sij may depend on Iq, reflecting 
evolution of the mass and spin. However, since the M and Sij in (16) do correspond to 



those defined in (11) and (13), we know that these quantities are in fact independent of 
time to, as shown in paper I and remarked above. Thus the near-zone background metric 
is identical (through order 0(l/f^)) at all to, which implies by the form of ^ that the 
near-zone perturbation is independent of t through 0(l/f^), becoming simply 

9^ = D,^{n) + H,^{n)\ + O + 1 O (i) . (17) 

Further simplification can now be made. The perturbation g^^^ must satisfy the (vacuum) 
linearized Einstein equation about the near-zone background g^^p. The quantity D^i, appears 
at order 0(l/f^) (and higher) in the linearized Einstein tensor, while H^y appears first at 



0{l/r^). However, at these orders only the first two terms in the background (16) will 



appear, so that the background is effectively Schwarzschild. Therefore, when considering 



only the explicitly displayed terms in the perturbation (17), we may use well known results 



T6t [T7] for perturbations of the Schwarzschild spacetime. In particular, since stationary 



perturbations scaling as f° and (as in ( 17)) are known to be £ = (spherically symmetric) 
up to gauge, the perturbation in ( [l7) ) is simply a change in mass, and a gauge may be chosen 
so that 

,^y^^v.o(^).fo(i), (18) 

where 5M is an arbitrary constant (which here may a priori depend on to), which I refer to 



as the perturbed mass of the body. An explicit formula may be given in analogy with (11 ), 



5M(to) = ^ lim / n'di^^f^n, (19) 
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with the caveat that this equation holds only in coordinates where D^j^ = C^jj, = 0. If one 
wishes to compute the mass in other coordinates (such as when the Lorenz gauge is used, 
and D^j/ is equal to the value of the "tail integral"), a more complicated expression (which 
must implicitly involve transforming to appropriate coordinates) must be derived. However, 
since the mass is an intrinsic property of a spacetime, there is no need to consider such 
coordinates in defining 6M and determining its evolution.^ Note that since the perturbation 



(17) is still stationary at 0(l/f ), it may be possible to define perturbed spin and deviation 
by an analogous procedure. However, these quantities appear at one order in A higher than 
pursued in this paper, and are not considered here. 

I now show that the perturbed mass does not evolve with time. With our previous 
coordinate choices, the second-order near-zone perturbation takes the from 

gf^ = B^i^jx'x^ + fE^, + K^, + ^2<5MV + ^(l/f) + i 0{l/f^) + ? 0{l/f^). (20) 
This perturbation satisfies the (vacuum) linearized Einstein equation about gf} (equation 



([T6|)) with effective sources constructed from (/^p (equation ([18|)). However, it is easy to 



see that these effective sources are 0(l/r^), tO{l/r^), and PO(l/r^), while the error in 
the linearized Einstein tensor introduced by including only the explicitly displayed terms in 



the second-order perturbation (10) is 0(l/r^), tO{l/r^), and PO(l/r^). Thus the effective 



source terms may be ignored, and denoting the linearized Einstein tensor of g'^^^ about g'^^p 
by QfiD, we have Qp^p = to the appropriate order, i.e., 

g-,, = o(i/f3) + to(i/f4) + Po(i/f5). (21) 

To determine the mass evolution 5M we focus on the i = part of Qpp. While the 
background cjp^p is not spherically symmetric due to the presence of the spin term, it is easy 
to see that to the relevant order in Qnp the spin term only contributes in product with the "B" 



term in the perturbation (20 ). However, since the far-zone background spacetime is assumed 
to be vacuum, its Riemann tensor may be decomposed into two rank-two symmetric, trace- 
free spatial tensors (the "electric" and "magnetic" parts, e.g. [1]), showing that the "B" 



A definition of perturbed mass was given in [19] , which appears to correspond to equation ( 19 1 applied 
in the Lorenz gauge. This definition would not be sensible within our framework. The mass defined in 
|19j was found to evolve with time. The conclusion that the perturbed mass evolves with time appears 
to be at odds with the analysis of [2D] 7 where it was found that energy conservation is satisfied under the 
assumption of no change in mass. 
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term is pure i = 2. Since the spin term is £ = 1, the combination can therefore make no 
contribution to the i = part of Qfip at the relevant order. Instead, the i = part is 
completely determined by the remaining terms in the perturbation, which "see" only the 
Schwarzschild metric. In particular only the i = parts of these terms may contribute to 
the i = part of Qfip, and we conclude that the £ = part of g^p must be a perturbation 



of the Schwarzschild spacetime to the relevant order. We may now use Zerilli's result [16] 
that i = perturbations of Schwarzschild simply shift the mass, being writable as a 1/f 
term plus a gauge transformation. Since 0(l/f) is higher than considered, we have that the 
£ = part of g^^p is pure gauge to the relevant order, 

(rE^, + K^, + -26M6^,) =Vi^ ^ + 0(l/f) + i 0(1/ f^) + ?0(l/f3), (22) 

M 

where V is the derivative operator compatible with Schwarzschild and is a vector field. 
We may now take a time derivative to find 

do Vi, ^ =Vi, = —V + 0(l/f2) + t 0(l/f3), (23) 

which immediately implies 6M = 0, since the mass perturbation 25M /fS^i^ is not pure 
gauge. 

Incorporating this result, we may now summarize the coordinate choices made in this 
section as 

mi 2M , M2 AnHuS^)i ^ f l\ 
g^ = + —5,. + ^{3v,. - t,U) ^ + ^ ( ^ ) (24) 



^(1) ^ — + o / ^ ] ^toi ^ ] (25) 



26M ^ _ / 1 



Jill/ — ^flU I ^ \ — 9 I ' \ — T 



9§ = B^^^x^x^ + fE,^ + K,^ + + t O (^1) + ? O , (26) 

with 5*01 = 0. By fixing the mass and spin terms to a standard form at all to and choosing 
the mass dipole to vanish for all Iq, the metric form has been made very simple, and all 
non-stationarity has been eliminated from the relevant orders in f. These properties make 
this gauge much simpler to use than the Lorenz gauge used in [5] and elsewhere. Rewritten 
in the far-zone, the perturbation series in these coordinates becomes 

h,. = gl'J = ^-^ + E,^r + Oir') (27) 



9^"^ = ^(3^.^ - ^ + —V + K,. + 0{r), (28) 

-,(0) 



with the background g^J still given by g. 
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IV. PARITY CONDITION 



The definition of center of mass adopted in section |TT] is based on the dipole moment of 
the time-time component of a stationary, asysmptotically flat metric. An alternative defi- 
nition of center of mass is given by equation (5.13) of Regge and Teitelboim [7], derived as 
the conserved quantity canonically conjugate to the asymptotic boost symmetry of asymp- 
totically flat general relativity.^ Like the Hamiltonian notion of mass (normally referred to 
as the ADM mass), this "Hamiltonian center of mass" involves only the spatial metric, and 
is more general in that it may be applied to time-dependent spacetimes in addition to the 
stationary spacetimes we consider. However, unlike the Hamiltonian notion of mass, the 
center of mass comes with an additional restriction: In order to ensure the existence of the 
integral defining the center of mass, Regge and Teitelboim impose a "parity condition" that 
the monopole (l/r) part of the spatial metric be even parity, Cij{n) = Cij{—n) in equation 
^7 This restriction is not necessary to define the mass dipole, which is finite for any metric 
of the form (Isl). 



Since the general metric form (16) satisfies the parity condition, we see that the parity 
condition is simply a coordinate condition in the context of stationary, asjTiiptotically flat 
spacetimes. Rotations and translations will automatically preserve the parity condition, 
while a supertranslation 5a;* = a*(n) + 0(l/f) must now satisfy a* = c* -|- S*(ri), with c* 
constant and odd parity, S*(— n) = — S*(n). This form provides a natural split between 
the "pure translation" part c* and a "pure supertranslation" part S*, which is odd-parity. It 



is easily checked from equation ( 15 ) that the mass dipole center of mass changes by d under 
this transformation, so that the parity condition removes its supertranslation dependence. 
In this case the transformation properties of the mass dipole center of mass agree with those 
established by Regge and Teitelboim for their center of mass; and since both notions give 



zero on the metric (16), the two definitions are equivalent in our context. Therefore, the 
question of which definition to use is simply the question of whether to impose the parity 
condition. 

One may take one of two alternative viewpoints on this matter. First, since the formal 



^ We note that a later formula due to Beig and O'Murchadha [T2] is equivalent given the parity condition. 
^ Regge and Teitelboim also impose a parity condition on the extrinsic curvature. However, this condition 
is not needed for the center of mass and plays no role in our analysis. 
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Hamiltonian analysis yields a center of mass formula that diverges (in general) in coordinates 
that violate the parity condition, one may argue that such coordinates are "too irregular" 
to admit a notion of center of mass, even if the mass dipole formula is finite. Alternatively, 
one may view the mass dipole formula as providing an extension of the Regge Teitelboim 
center of mass to a larger class of coordinates within the stationary case. In any case, the 
parity condition adds a number of simplifying properties in the context of the present work: 
It eliminates the supertranslation dependence of the mass dipole, it allows the equations 
of motion to be expressed purely in terms of the local spacetime metric (see discussion in 
appendix [B]) , and it makes the mode-sum regularization scheme gauge invariant. 



V. EQUATION OF MOTION 



In the gauge of section III, the equation of motion for the deviation is simply Z^{t) = 0. In 
principle, therefore, giving the change in under a change in gauge provides the complete 
description of motion. However, the more useful description of motion in other gauges, 
equation ([T]), may be derived as follows. Begin with gauge transformations x'^ = + X^^ 
of the form® 

e = a^(t,n) + 0(r), (29) 

where ^'^ is assumed smooth in r at fixed {t, n), so that in particular a'^ is a smooth function 
of its arguments. It is easy to check that such transformations preserve the form of equations 
(|5]|6]), and thus are allowed under our assumptions. As described in [T3], we furthermore 



believe (but have not proven) that expressibility in the form of equation (29) is a necessary 
condition on an allowed transformation, except in the case of certain trivial one-parameter- 
families, where additional logr terms are allowed. Thus we believe (but have not proven) 
that such transformations correspond precisely to the coordinate choices allowed by our 
formalism (not including the parity condition) at first order in A, for non-trivial families of 
solutions. In order to preserve the parity condition Cij(t, —n) = Cij{t,n), we must restrict 
the form of so that 

e = c\t) + J:\t,n) + Oir), (30) 

® In the appendix of paper I an opposite sign convention, x'^ — x'^ ~ X^^ , was used for the definition of the 
gauge vector. 
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with odd-parity, S*(t, — n) = — S*(t,n). I will refer to gauge transformations of the 
above form as parity- regular transformations, and I will define parity-regular gauges as 



those that are related to the gauge of section III by a parity-regular transformation. We 
believe that parity-regular gauges are the general class allowed by our assumptions plus the 
parity condition, except possibly in trivial cases. Thus we believe that one may check if a 
given gauge is parity-regular by checking that the metric perturbation takes the form ^ 
with Cij{t, —n) = Cij{t,n). However, absent a complete proof of the assertions in [H], one 
must instead check that the gauge vector relating to a known-parity-regular gauge is of the 



form ( 30 ) . At the end of this section I discuss the parity-regularity of some common gauge 
choices. 



Under a change of gauge of the form ( 30 ) , the near-zone background coordinates change 



by x'* = + d{to) + S*(to,?^) + 0(l/f). Using equation (15), we see that the deviation 
Z^lto) changes by c*(to), as noted in the previous section. We may express this in terms of 
the gauge vector by taking an angle-average over a small constant-r sphere,^ 

SZ' = iOr^o = ^ lim f Cdn. (31) 

47r r-¥0 / 



Equation (31 ) gives the change in deviation due to a parity-regular transformation made on 
any perturbation of the assumed form ([6]). 

The key manipulation now follows. Consider the second time derivative, 6Z^ = dodoSZ^. 
We have 

6Zi = (dodo^,) 

= (VoVo^^ + Rojokx'^dj^i + VoVi^o - ViVo^o + Rioo'Qr-^o 

= 6{Fi)r^o - Rm.^Z^ (32) 

In the second line we have rewritten in terms of covariant derivatives {Romx'' dk^t is a 
Christoffel term), as well as added zero in the form of the Ricci identity. However, since 
dj^i is even parity to leading order and the Riemann tensor is smooth, the first Riemann 



^ We could equivalently express 5Z^ as an average over a circle or over two antipodal points, since these 
averages all agree for a "constant plus odd parity" function. The entire derivation of the equation of 
motion could then proceed unchanged, so that the self-force in fact may equivalently be written as the 
average of the bare force over a (constant geodesic distance) sphere, circle, or pair of points. 
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term vanishes by virtue of the parity condition. Noting that the remaining derivatives of 
appear only in symmetrized form, in the third hne we reexpress in terms of the change 
in the metric perturbation, = — 2V(^^i,), finding precisely the gravitational force form 
of equation (|2|. We also separate off the remaining Riemann term, which takes a geodesic 
deviation form. In the last line we use equations ^ and (31), where the finiteness^° of 



(Fj)r-s.o allows us to pull the 6 out of the angle- average. We may now "drop the 5's" to 
obtain 

Z' = {F')r-,o-Rojo'Z'+A\ (33) 

where is the constant of integration — an unknown gauge-invariant acceleration. Thus 
by simple manipulations we have fixed the form of the equation of motion up to a gauge- 
invariant piece, and may now work in any convenient (parity-regular) gauge to determine 



Since the gauge of section III has Z'''(t) = it immediately eliminates two terms in (33), 
giving simply 

A' = -(F%^o (34) 

in this gauge. The interpretation of the gauge is that the gravitational self-force is exactly 
opposite to the gauge-invariant force, so that the total force is zero, and there is no deviation 
from geodesic motion {Z'^(t) = 0). Since the angle-average of a three- vector picks out an 
i = 1, electric parity part, we need consider only the i = 1, electric parity part of the 



bare force in order to compute A^ from equation (34). Since only the "E" term in the 



perturbation (27) contributes, we may focus on the i = 1, electric parity part of E^^. To do 
so, we return to second-order near-zone perturbation theory. 

As remarked above in the derivation of the constancy of the perturbed mass, the relevant 

_(2) 

terms in the second-order metric perturbation (/^p satisfy the linearized Einstein equation 



off of the gfp to the relevant order, equation (23). As further remarked, the £ = 1 spin 



term in the background appears only in combination with the £ = 2 "B" term in the 
perturbation. While there is no contribution to the £ = mode relevant for the perturbed 
mass, there can be a contribution to the £ = 1 mode relevant for the deviation, which 
significantly complicates the analysis (see appendix IaI) . However, if the spin is assumed to 



10 



The angle average of is manifestly finite in the gauge of section III and, by reversing the calculations 



of (32), may be easily seen to transform finitely (see also (40l) 
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be zero from the outset, then, as in the mass evolution case, the metric perturbation "sees" 
only Schwarzschild to the relevant order, and we can make use of Zerilli's analysis of 
perturbations of Schwarzschild. In particular, Zerilli showed that i = 1, electric parity 
perturbations are pure gauge, so that by a (second-order near-zone) gauge transformation 
we may eliminate the i = 1, electric parity part of Ef^^, entirely, whence it immediately 



follows from (34) that = 0. Thus the equation of motion in the spinless case may be 



derived with very little effort, involving only the few lines of algebra of equation (32). If 
the spin is not assumed zero, more algebra is required (though still significantly less than 
needed when the Lorenz gauge is used) . This case is treated in appendix |A| leading to 

MA, = ^S'^Rkm, (35) 

showing that A^ is simply the acceleration due to the Papapetrou spin force. We have thus 
justified the final equation of motion, which appears in covariant form in ([T]). 

I now discuss the parity-regularity of some common gauge choices. The Lorenz gauge is 



convenient for local series expansions about the particle. From equation (|27|), the gauge of 
section 



III already satisfies the Lorenz condition at leading and subleading order, V^{h^y — 
{l/2)hg^^) = 0(1). It is then easy to check that the gauge vector to a full Lorenz gauge 
must take the form ^' = C*(t) + 0(r), so that the Lorenz gauge is parity-regular. 

The Regge- Wheeler gauge is convenient for perturbations of Schwarzschild. Barack and 



Ori [6] show that in a few specific cases the gauge vector takes the form (29) (bounded but 



direction dependent), and their formulae also imply that in these cases one can choose the 



vector to satisfy the required parity property, equation (30). However, Hopper and Evans 
[TH] have shown that in general the Regge- Wheeler gauge metric perturbations contain a 
delta-function in the Schwarzschild radial coordinate at the position of the particle, so that 
in general the Regge- Wheeler gauge is too singular to define the motion by our procedure. 
However, using the explicit formulae for the coefficient of the delta function found in [15], it 
should be possible to simply eliminate the delta function by a gauge transformation during 
the process of reconstructing the metric perturbations. If the resulting gauge is parity- 
regular (a suggestion consistent with Barack and Ori's restricted results), then the results 
of this paper would enable self-force computations to be made within the Regge- Wheeler 
formalism, using mode sum regularization if desired. 

The radiation gauge is convenient for perturbations of Kerr. In [9] it was checked that 
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a radiation gauge may be chosen so that the metric perturbations near the particle take 
the form {l/r)C^^{t,n), with Cij{t,n) even parity on the sphere, so that the singularity 
is properly 1/r and the parity condition is satisfied. However, it was not checked that 
Ctiu{t,n) is smooth (as assumed in this paper), and in fact it can be seen from the analysis 
of [in] that C^iy(t,n) contains a jump discontinuity along a great circle, which is inherited 
from a discontinuity in the metric perturbation located at the radial coordinate of the 
particle. Such a discontinuity seems unlikely to threaten the validity of the results, since all 
of the relevant formulae remain defined. However, in order to be certain that the results of 
this paper may be applied to the gauge of a careful analysis of the gauge vector relating 
to some parity-regular gauge must be performed. Armed with the explicit form of such a 
vector, it should be straightforward to check if the arguments of this paper still hold. 



VI. MODE-SUM REGULARIZATION 



The computation of gravitational self-forces on black hole background spacetimes is an 
important problem for gravitational- wave astronomy of extreme mass-ratio inspirals (e.g., 
|21j). The angle- average formula suggests a straightforward way of proceeding: first numer- 
ically compute the metric perturbations of a point particle in any parity-regular gauge, and 
then perform an average to determine the force in that gauge. However, while simple in prin- 
ciple, such a procedure may be difficult to carry out in practice, due to the singular nature 
of the quantity being averaged. Instead, to achieve an accurate result it is likely preferable 
to use an alternative technique, such as that of mode sum regularization, first introduced in 
[8] and widely employed thereafter. This method takes advantage of the fact that numerical 
calculations in black hole spacetimes usually employ a spherical (or spheroidal) harmonic 
decomposition, which in particular has the property that the individual modes of the bare 
force are finite at the particle. One then performs a finite subtraction on each mode, which 
is designed so that the resulting sum converges to the correct self-force. Extensive work 
has determined the form of this subtraction (in terms of "regularization parameters") for 
arbitrary orbits of Schwarzschild and Kerr in the Lorenz gauge. In this section I show that 
the mode sum regularization scheme is gauge-invariant under the parity condition, in the 

At the time of the writing of [9 , I beUeved (and communicated to the authors of [5]) that only this check 
was required for the results of the present paper to hold. 
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sense that the same subtraction may be employed to determine the force in any parity- 
regular gauge. Combined with the Lorenz gauge results of [221, ^^^^ provides a complete 
regularization prescription for Kerr in parity-regular gauges. 

Let (i, f, 9, 0) be Boyer-Lindquist coordinates for the Kerr spacetime. For a given point 
along the worldline 7 where we wish to compute the self-force, we may rotate and time- 
translate the coordinates so that the particle is located at t = = 0, while taking the 
remaining coordinate positions to be f = tq and 9 = 9q. Despite the lack of a full rota- 
tional symmetry we (following [E]) nevertheless also perform an additional rotation in the 9 
direction, so that the particle is located at the pole of the new coordinates. More precisely, 
define new angular coordinates {9, 0) by 

cos ^ = — sin 6' cos cos 9o + cos 9 cos 9o (36) 

tan = , , (37) 

sin 9 cos cos 6*0 + cos 9 sin 9o 

to obtain "rotated Boyer-Lindquist" coordinates t, r, 9, in which the particle position is 
given by t = 0, f = To, 6' = 0. In these coordinates the metric components are smooth every- 
where except for the pole 9 = 0, where they acquire non-trivial direction-dependent limits. 
Below we will need the lowest-order relationship between the spatial Fermi coordinates 
and the rotated Boyer-Lindquist coordinates, restricted to the sphere. A straightforward 
computation gives this to be 

x^|t=o,f=ro = a* ^cos0 + /3' 9 sin 4) + 0{9^), (38) 

where a* and are constants independent of 9 and (dependent on the Boyer-Lindquist 
position the three- velocity, and the mass and spin parameters of the Kerr metric), and 
where the 0(9'^) term may depend on 0. 

The advantage of placing the particle at the pole is the simplification of the spherical 
harmonic description by the elimination of modes with non-zero m when the series is evalu- 
ated at the particle. Let a subscript i denote the tth term in the expansion evaluated at the 
pole/particle, = J AY^QdCl for integrable functions A{9,(f)), and suitably generalized for 
distributions. When viewed in light of our angle-average result, the mode sum regularization 

prescription amounts to finding some such that 

00 

(nr^o = i^i - Si) ■ (39) 

17 




FIG. 1: A diagram illustrating the geometrical setup of the mode sum regularization argument. 
The particle is at the pole of rotated Boyer-Lindquist coordinates. The mode decomposition is 
taken relative to the background coordinate sphere, while the self-force is given by an average over 
the local inertial sphere. For the change in bare force, a general theorem relates the mode sum at 
the particle/pole to the average over the polar circle, which agrees with the average over the local 
inertial sphere when the parity condition is satisfied. 

This formula relates an average over an infinitesimal sphere surrounding the particle to a 
spherical harmonic decomposition on a finite sphere surrounding the black hole, evaluated 
at the particle (see figure [T]). Such a connection between mode sums and local averaging is 
familiar from ordinary Fourier series, where, if a function is of bounded variation, its series 
converges to the two-sided average at a discontinuity. For spherical harmonic expansions, 
an analogous result (e.g. section 11122b of [23J) states that if the average of a function over 
latitude lines is of bounded variation (as a function of latitude), then its spherical harmonic 
series evaluated at the pole ("Laplace series") converges to the average on an infinitesimal 
latitude line surrounding the pole. (This result is easily understood at a formal level, by 
noting that Y^o is independent of 0, so that the formula for Ag takes an average over 0.) 
The theorem does not apply to the bare force (which is divergent), but it does apply to 
the change in bare force under a change in gauge, 6F\ In particular, a simple calculation 



(which reverses the calculations internal to the average in (32)) gives 



6F^ = dodoC + Rooj^' - RoVdkC + 0(r) (40) 



for any transformation of the form (29). If the transformation is parity-regular, equation 
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(30), we see that has a Fermi coordinate expansion of the form 



SF' = C{t)+S\t,n) + Oir) 



(41) 



where 5* is smooth and odd-parity, S'^{t,—n) = —S'^{t,n)}'^ Restricting to the background 
coordinate sphere t = 0,r = tq and expanding in 6 at fixed 0, we have 



\i=0,f=fo^ 



C^(0) +5^(O,n|,-=o,,=,o(^~,0)) + O(0~) 
C\0) + S\0, n\i=o,f=ro (0 = 0, 0)) + 0(6) 



where the 0{9) terms may depend on cp and we have defined 



no{(j)) = \imn(t = 0, r = tq, 9, 

9^0 



(42) 
(43) 
(44) 

(45) 



In moving from equation (42) to equation (43) we have used the fact that the restriction 



of n to the background coordinate sphere is smooth in 6 at fixed (see equation (38) and 



recall n* = x^/r), as well as the fact that 5* is smooth in n. Equation (|44j) shows that the 
restriction of to the background coordinate sphere is continuous (in fact smooth, by our 
assumptions) in 9 at fixed (p. In particular its average over (p is of bounded variation, and 
by the theorem we have 

1 



lim 



1 

2^ 



t=0,f=rol^5 



(46) 
(47) 



where in the second line we have plugged in the form of equation (44). Since 5* is odd 
parity, the second term on the right-hand-side will vanish if no is odd under — )■ + 7r. This 
property is expected from the geometry of the setup (figure [T|, and is easily confirmed from 



equations (45) and (38). Thus the term involving iS* vanishes, so that the Laplace series for 



SF^ in fact converges to C*. However, the angle-average that computes the self-force also 



returns C* on the form (|41|). Therefore, when the parity condition is satisfied the averages 

J2iSF% = {5F%^o, (48) 



agree, and we have simply 



The parity condition is not required to show the apphcabihty of the theorem, but will be necessary for 
the later analysis of this section. 
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showing that the Laplace series for the change in bare force 5F'^ in fact converges to its local 
inertial angle- average, i.e., to the change in self- force it effects. This means in particular 
that no extra must be subtracted in the new gauge, since merely the process of decom- 
posing into modes and resumming returns its contribution to the self-force. To see this 
explicitly, let F^^^ denote the bare force in a gauge that satisfies the parity condition, and 
write 

= + 6F^)r^,. (49) 



In writing the second line we have used (48), and in writing the third line we have used 



equation (39). This shows that equation (39) holds in the new gauge if it held in the old, 



i.e., that S\ is a correct piece to subtract in any parity- regular gauge. 

We note that previous work has organized the subtraction so that one first subtracts an 
5"^ of the form S\ = A\l + 1/2) + 5* + C'/{1 +1/2), then sums over modes (the result is 
now finite), and finally adds in a finite residual to get the correct self-force. (In terms of 
our S\, is a "finite piece" = J^eiS} - S})}^) The (^-independent) A,B,C,D are the 
"regularization parameters" for the particular orbit and spacetime, which by the results of 
this section do not depend on the choice of parity-regular gauge. 
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Remarkably, it has been found (by lengthy computation in the Lorenz gauge) that Z?* = in every 
circumstance, so that the subtraction of A'-{£ + 1/2) + _B' + C^/{£ + 1/2) in fact returns the correct 
force. This surprising relationship between the large-i expansion of a point particle metric perturbation 
(which uniquely determines A, B, C) and the physical self-force has thus far defied a more fundamental 
explanation. 
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Appendix A: Derivation in the case of non-zero spin 



When the spin is non-zero, we may not rely on Zerilh's results for Schwarzschild [TB] , and 



the analysis of second-order Einstein equation, equation (23), becomes more complicated 



In this case it pays to systematically consider the contributions to the linearized Einstein 



tensor Qpp from the various terms in the background (24) and perturbation (26). Since 
all terms are stationary to the relevant orders, no t-dependence will appear, and we may 
count orders in 1/f. At leading order 0(1) in Q^p, the only contributions are from the "B" 
term in the perturbation ( [26] ) and the flat "r/" term in the background. Since the "B" term 
is (by the Fermi coordinate construction) a perturbation of fiat spacetime, the linearized 
Einstein equation is automatically satisfied and we learn no new information. At next order 
0(l/f) in Qfip, both the "E" term in the perturbation and the mass term 2M5^p/f in the 
background can now contribute. Expanding the background metric in powers of 1/f, the 
linearized Einstein equation at 0(l/f) may be written 

Gji) [E^^f] + 2Gf [2M/f5^„ B^,,^x'x^] = 0, (Al) 

where G^^^ and G^^^ are the first and second-order Einstein tensors (respectively) off of flat 
spacetime. However, since the "B" term has no £ = 1, electric parity part while the M 



term is spherically symmetric, the second term in (Al) has no £ = 1, electric parity part. 
Therefore the £ = 1, electric parity part of the "E" term satisfies the vacuum linearized 
Einstein equation about fiat spacetime, 

GW[(i?,,,f),=i,+]=0, (A2) 

where the subscript = 1, -|-" indicates the i = I, electric parity part. Since i = 1, electric 
parity perturbations of fiat spacetime that scale linearly with f are pure gauge (e.g., [IS]), all 



solutions to equation (A2) may be written (_E'^j,f)£=i + = d(^^£y), where £^is ^ = 1 and electric 
parity. If we make a (second-order near- zone) gauge transformation generated by then we 
may set the ^ = 1, electric parity part of E^^ to zero. However, one can check that in general 
the required gauge vector is t-dependent (despite the "E" term being stationary), so that 
this gauge transformation would introduce t-dependence at higher orders, in violation of our 
previous choices that eliminated such dependence. (In particular, terms of order i/f and 
P/f^ would appear in gp^^, contradicting the previous choices H^^i, = C^y = — see equations 
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(10) and (26).) Without invalidating these choices we may only make gauge transformations 
whose gauge vector is t-independent. One may check that this allows us to put E^^, in the 
form 

{E^^f)i=i^+ = -2aixH^U, (A3) 

where a, is an arbitrary spatial vector, named since in this form {E^yf)i=i^j^ is an "accelera- 
tion perturbation" familiar from Fermi coordinates about an accelerated worldline. Making 
such a gauge transformation (and "absorbing" its effects at 0(1) into the arbitrary tensor 



K^i,), equation (34) becomes simply 



A' 



a . 



(A4) 



We have now made coordinate choices that eliminate all t-dependence and reduce the 
relevant ^ = 1, electric parity part of E^^, into a simple form with one unknown, the gauge- 
invariant acceleration = . It remains to use the linearized Einstein equation at order 
0(l/f^) to determine A^. Again expanding the background in powers of 1/f, at this order 
the = 1, electric parity parts of the second-order Einstein equation may be written as 



[K^,] + 2Gf ) 



2M5. 



+ 2Gi2) 



B 



1 ^fj.iiyj 



0, 

(A5) 



where terms that can give no i = 1, electric part have not been displayed. This equation 
gives relationships between Kfj,^, Ai, M, Sij, and Rf^iuj. To determine a relationship not 
involving the unknown tensor K^^, first write out the linearized Einstein tensor about flat 
spacetime for stationary perturbations, and note that there is a particular i = 1 part^^ 
that vanishes for all K^,^. Then computing this particular i = 1 part of the remaining two 
terms will determine Ai in terms of the mass, spin, and Riemann tensor. Performing this 
calculation yields MAi 



^S^'-Rkioi, as claimed in equation (35). 



Appendix B: Equation of motion in parity-irregular gauges 



I now consider the form of the equations of motion in parity-irregular gauges, adopting 



the mass dipole definition of center of mass. Under the general gauge transformations (29), 



In the notation of [5], this part is {2/3)Rf + Rf, involving only spatial trace-free components of the 
linearized Ricci tensor. 
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we have from equation ( 15 ) that^^ 

5Z' = — lim / i^nm' = 3(f n.n^^o- (Bl) 
47r r-s>o J 

Beginning with the equation of motion in a parity-regular gauge, we may now derive an 
equation for in a new gauge, 

{{-dodo^d^, - 3n\) + Rojokx'd'C - RJU^\ - 3n'=n,))),,^o, (B2) 



where equation (40 ) has been used. If the parity condition is satisfied, and ^* = c* + + 
0(r) with odd-parity, then the right hand side vanishes and the equation of motion retains 
the original form, depending only on the local spacetime metric (at zeroth and first order in 
perturbation theory). However, if the parity condition is not satisfied, then the right hand 
side does not in general vanish, and the equation of motion for Z'^ takes a complicated 
form involving the gauge transformation to some reference gauge. Another way to see the 



difficulty is to repeat the calculations of (32) for a general gauge transformation, giving 

Without the parity condition the Riemann terms do not simplify into the geodesic deviation 
form Rf)iQj5ZK It appears that no expression in terms of just and dh^y is possible, so 
that the equation for in parity-irregular gauges must involve a gauge vector explicitly. In 
particular, there appears to be no natural separation of the terms contributing to into 
"self" and other forces in the case of a parity-irregular gauge. 
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